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Abstract: We study the thermoelectric conductivities of a strongly correlated system 
in the presence of a magnetic field by the gauge/gravity duality. We consider a class of 
Einstein-Maxwell-Dilaton theories with axion fields imposing momentum relaxation. Gen¬ 
eral analytic formulas for the direct current (DC) conductivities and the Nernst signal are 
derived in terms of the black hole horizon data. For an explicit model study, we anal¬ 
yse in detail the dyonic black hole modified by momentum relaxation. In this model, 
for small momentum relaxation, the Nernst signal shows a bell-shaped dependence on the 
magnetic field, which is a feature of the normal phase of cuprates. We compute all alternat¬ 
ing current (AC) electric, thermoelectric, and thermal conductivities by numerical analysis 
and confirm that their zero frequency limits precisely reproduce our analytic DC formu¬ 
las, which is a non-trivial consistency check of our methods. We discuss the momentum 
relaxation effects on the conductivities including cyclotron resonance poles. 
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1 Introduction 

Strongly coupled electron systems show many interesting phases such as non-Fermi liquid, 
high Tc superconductor and pseudo gap phase. Some of the most important and basic 
experimental observables in investigating those systems are the conductivities: electric((T), 
thermoelectric(a, d), and thermal(k) conductivity. Therefore, it is essential to develop a 
theoretical method to compute conductivities to explain and guide experiments. However, 
due to strong coupling, the perturbative analysis of quantum field theory does not work 
and we don’t have a reliable systematic tool to compute them. 

Gauge/gravity duality is an approach for such strong coupling problems, and it has 
been developed as a method for conductivity [1, 2]. Some early works treated systems 
which have translation invariance. However, at finite charge density the direct current (DC) 
conductivities in such systems are infinite. To solve this problem, it is essential to introduce 
momentum relaxation. For this, several ideas have been proposed. 

The most straightforward way is to impose inhomogeneous boundary conditions of the 
bulk fields to break translation invariance [3-8]. Massive gravity models studied in [9- 
13] give mass terms to gravitons, which break the spatial diffeomorphisms (not the radial 
and temporal ones), and consequently break momentum conservation of the boundary 
field theory with translational invariance unbroken [11]. Holographic Q-lattice models and 
those with massless linear dilaton/axion fields [14-23] take the advantage of a continuous 
global symmetry of the bulk theory. Some other models utilise a Bianchi VHq symmetry 
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to construct black holes dual to helical lattices [24-27]. All these models yield finite DC 
conductivities as desired. 

However, all models with momentum relaxation, except [23], did not include the mag¬ 
netic field. Since the transport properties at finite magnetic field such as the quantum Hall 
effect, the Nernst effect, and the Hall angle, are also important basic probes for strongly 
correlated electron system, it is timely and essential to develop the methods for them in the 
presence of momentum relaxation. Indeed, the holographic analysis on conductivities at fi¬ 
nite magnetic held was one of the pioneering themes opening up the AdS/CMT (condensed 
matter theory) [28-30]. The purpose of our paper is to extend them by implementing mo- 
menturm relaxation holographically^. This paper is also a companion of [19, 32, 33], where 
all DC/AC electric(cj), thermoelectric(Q:,d), and thermal(k) conductivities are analysed 
thoroughly in the absence of magnetic held. 

We consider a general class of Einstein-Maxwell-dilaton theories with axion helds im¬ 
posing momentum relaxation. First, we derive the analytic formulas for DC electric((T), 
thermoelectric(a, d), and thermal(k) conductivity in terms of black hole horizon data fol¬ 
lowing the method developed in [34]. Based on these formulas we discuss the model inde¬ 
pendent features of the Nernst signal. Notice that the Nernst signal (2.55) 

e-N = —(o' • a)y , 

is zero in the holographic model without momentum relaxation, since the electric conduc¬ 
tivity is inhnite. Thus, momentum relaxation is essential for the Nernst effect. The Nernst 
signal has interesting properties which could support the existence of the quantum critical 
point(QCP). As we approach to the QCP or the superconducting domain, the strength of 
the Nernst signal becomes stronger and shows a non-linear dependence on the magnetic 
field, which is different from the expectation based on the Fermi liquid theory [35] See 
[29, 31], for pioneering works on the Nernst effect by the holographic approach and the 
magnetohydrodynamics with a small impurity effect. We deal with similar topics by means 
of a general class of holographic models encoding momentum relaxation, where we assume 
that momentum relaxation is related to finite impurity density, which could be large. Note, 
however, that this relation is not proven yet. 

After discussions on a class of models, we study in detail the dyonic black hole back¬ 
ground [28-30], modified by the specific axion fields introduced in [16]. We numerically 
compute AC electric(iT), thermoelectric(Q!, d), and thermal(k) conductivity and confirm 
their zero frequency limits agree to the DC formulas that we have derived analytically. It 
recovers the results in [28] if momentum relaxation vanishes. We discuss the momentum 
relaxation effect on the conductivities including the cyclotron resonance poles, which was 
first observed in [30]. 

This paper is organized as follows. In section 2, we consider a general class of Finstein- 
Maxwell-dilaton theories with axion fields and derive general formulas for the DC electric, 

^In [31], momentum relaxation is introduced perturbatively in the hydrodynamic limit. 

^These anomalous behavior may be described by a liquid of quantized vortices and anti-vortices in 
non-superconducting phase. See [37] for a speculative point of view. 
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thermoelectric, and thermal conductivity at finite magnetic field as well as the Nernst 
signal. In section 3, as an explicit example, we analyse the dyonic black brane with the axion 
hair and discuss the Hall angle and the Nernst effect. In section 4, we continue our analysis 
on the model introduced in section 3. We compute the AC electric, thermoelectric, and 
thermal conductivity numerically. The momentum relaxation effect on AC conductivities 
and the cyclotron resonance poles are discussed. We compare the zero frequency limit of 
our numerical AC conductivities with the DC analytic formulas derived in section 3. In 
section 5 we conclude. 

Note added: While this work was near completion, we noticed the appearance of 
[38-40] which have some overlap with ours. [38] deals with a massive gravity model at 
finite magnetic field. [39] considers the same class of models as ours. [40] obtains general 
expressions for conductivities at finite magnetic field using the memory matrix formalism. 


2 General analytic DC conductivities at finite magnetic field 


In this section, we will derive analytic formulas for the DC conductivities((T, a, a, k) in the 
presence of a magnetic field, from a general class of Einstein-Maxwell-Dilaton theories with 
axion fields(yi, X 2 ) 


S= 


R-1, [{d4>f + + ^2mdX2?] - V{ct>) - 






( 2 . 1 ) 


where F = dA and all 4>i ((/>), <^2 (</>), and Z{(l)) should be non-negative for the positive en¬ 
ergy condition. Without the magnetic field, analytic formulas for conductivities((T, a, d, R) 
were computed in [34] while, with the magnetic field, the electric conductivity a was pre¬ 
sented in [23]. Here we compute all the other conductivities as well. We employ the 
method developed in [8, 34], but a finite magnetic field poses some technical subtlety. We 
will explain how to treat it. The action (2.1) yields equations of motion: 


RmN - - x(dM'/>)(dAr(/>) “ X ^ ^i{4‘){dMXi){dNXi) - ^^RuRpN =0 , (2.2) 


i=l,2 




Vm (^(-/>)T^'^) =0 , (2.3) 
VM(‘l'i(<^)V“xi) =0, (2.4) 

dZ{4>) 


dcj) 


=0, (2.5) 


where, £ is the Lagrangian density of the action (2.1). 

To study the system at finite chemical potential with a background magnetic field(H), 
we take the gauge potential as 

A = Am^lx^ = a{r)dt + —{xdy — ydx). (2-6) 
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We choose the axion helds 


Xi = hx , X2 = k2y , 


(2.7) 


which break translational invariance and can give rise to momentum relaxation [14, 15]. A 
metric anasatz consistent with the choice (2.6) and (2.7) is 

ds2 = GuNdx^dx^ = -U{r)dt^ + j^dr'^ + e^dr)^x^ + ^ ^ 2 . 8 ) 

U[r) 

with (j) = 0(r). We consider the case in which the black hole solution exists and its horizon 
is located at r = r^, i.e. U{rh) = 0. We further assume the ansatz 

4>i(</)) = $ 2 (</>)=: ^((/>), vi{r) = V 2 {r) =■. v{r), ki = k 2 =: P . (2.9) 

First, the Maxwell equation (2.3) yields a conserved charge p in the r-direction 

p = = Z(0)e"Ma'(r), (2.10) 

which is identihed with the number density in the boundary held theory. The axion equa¬ 
tion (2.4) is trivially satished. The Einstein equation (2.2) and scalar equation (2.5) become 


v'\r) = - {rf + (l)\rf) , 


( 2 . 11 ) 


U'{r)=- 


^-2v{r) 


2Z{(l))v'{r 
_l_ e^(’')[/(r)(n'(r)^ — 4>'{r)‘^) 

-2v(r).2 


p2 ^ 2e2^WE((^)Z(0) + + e^^'''>Z{(l)){2/3‘^^{(j)) 


( 2 . 12 ) 


^—2v[r)^2 TT(r\ 

U"{t) + i?'e-2-Wz(</.) + e-"W/32cD(</,) + ^ {v\rf - , (2.13) 

by which we can obtain the background solutions for given Z{(f)) and V{(f)). For example, 
if we choose 


4>{r) = constant. 


ck( 0 ) = 1 , Zi^) = 1 


V(c/>) = -6/L^ 


(2.14) 


the background becomes the AdS-dynonic black hole geometry with the momentum relax¬ 
ation. We will discuss it in section 3 in detail. 

To compute the conductivities for the general background, we consider small pertur¬ 
bations around the background obtained by (2.11)-(2.13) 


6Ai = t6fl^\r) + 6ai(r) , 

(2.15) 

SGu = t5fP{r) + 5gti{r ), 

(2.16) 

SGri = e"(^)5(7«(r), 

(2.17) 

^Xi = ^Xi(r), 

(2.18) 
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(2.19) 

( 2 . 20 ) 


where i = 1,2 and t denotes the time coordinate. Sf^^\r) and are chosen 

= -Ei + Cia{r), 

Sf\r) = -QU{r), 

where Ei is the electric field and —Q is identified with the temperature gradient [34]. In 
spite of the explicit t dependence in (2.15) and (2.16) all equations of motion of fluctuations 
turn out to be time-independent, which is the reason to introduce the specific forms of 
(2.19) and (2.20). Furthermore, the electric current and heat current can be computed as 
the boundary(r —>■ oo) values of J*(r) and Q^{r), 

J\r) = Z{^)^E^^ir) 

= -Z{4>) {U{r){-e^^B5grj{r) + 5^^ 5a'^{r)) + a'{r)E^ 5gtj] , 

Q\t) = U\r) - ^ir)J\r ), 

which were identified in [34] at i? = 0 and they are still valid at finite B. 

Our task is to plug the solutions of the fluctuation equations into (2.21) and (2.23) 
and read off the coefficients of Ei and Q. First, the Maxwell equation yields 

0 = Om {Z{ct>)^gr^) 

= dr {Zi<P)^E^^) + dt {Z{cP)^r^) (2.24) 

= drB - Bdh-^^^\jZ{(t)). 

Therefore, the current at the boundary is given by 

poo 

/(oo) = /(r/ + Bd^Cj / dr'e-"("')z(,/.(r')) 

Jrh 

= f{rh) + Bd^CjEi. (2.25) 

Next, let us turn to the heat current, Q*, (2.23). It is convenient to start with the derivative 
of Q\ 

drQ^ = 2U{r)S^^dr - «'(// - «(//' 

= - 2Bd^CAr)e-<^^Z{(f). (2.26) 


( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 
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After using the Einstein equations for fluctuations with the ansatz (2.17)-(2.18), 


2U{r)5gl{r) = [2B ^+ 2/32ch(<i6)e-"M + U{r){v'\r) - (t)'\r))]5gt^{r) 

- 2pe-"(’')t/(r) {6a'^{r) + B5gry{r)] + 2BZ{^)e-<''^ {Cyair) - Ey} , 
2U{r)6g'ly{r) = {2S2z(</.)e-2"W + 2/32$(,/,)e-"M + U{r){v'\r) - ct>'\r))^ 5gty{r) 

- 2pe-<^'>U{r) {6a'y{r) - B5gr:c{r)} - 2BZ{^)e-<^^ {Cxa(r) - E^} , 


^Qrx — 


Bpe-<^^6gty + BZ{ct>)U{r)5a'{r) 


t/(r) (S2Z((/)) + ^2gKr)$(^)) L '--y' 

+ k^{cj))U{r)e^^'''^6xx{f') — pEx + U{r)v\r)'j Ci 

= C/(r)(B=ZW)+,S2e-M$W) [ “ - BZm(r)SK(r) 

+ k<^>{(j))U{r)e^^^^5x'ir) - pEy + (pa{r) - e^W[/'(r) + U{r)v'(r)) Cy 


(2.27) 


we end up with a relatively simple expression for the heat current; 


g*(oo) = Q\rh) + Be^^Ej / dr'Z{(t){r')) - 2Be'Kj / dr'Z{(l){r')) 


/ro 


'ro 


= Q*(r;,) + Be'^EjEi + • 


(2.28) 


In summary, we have two boundary currents: 


J'{oo) = J'{rh) + B^^CjEi, (2.29) 

Q'ioo) = Q\rh) + Be'^EjEi + , (2.30) 


where Q'irh) and J'irh) are functions at horizon, which can be further simplified by the 
regularity condition at the black hole horizon [8, 34] 


5ai{r) -ln(r - r/,) + • • • , 

5gti{r) ~ + 0{{r - rh)) ^ -, 


^Xi{r) ~ + 0{{r - rh)) + 


(2.31) 


Thus, the boundary currents yield 

J\oo) = -{pd'^e-"'^ + Be'h-'"^Zh)6g[f + 5'^EjZh + , (2.32) 

Q\oo) = -A-KT5'^5g^^'^ + Be'^EjJ^i + ^e'^Ci^a , (2.33) 

where Vh = v{rh) and Z^ = Z{(j){rfi)). 5g\^^ also can be replaced by the horizon data using 
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the equations of motion (2.27). The near horizon expansion of the last two equations in 
(2.27) gives 


X (^) 

^9tx = 

X (^) 

^9ty = 


o^h 


B^Zh + 

(,Vh 




Bpe-^'^6giy^ - BZ^Ey - pE^ - 47re"'*rC. 


- Bpe-^^5gt^ - BZhE^ - pEy - 


which, in turn, give us two algebraic equations for 6g^^ of which solutions are 


^9u = 


pt’h 


B2^p2 + 52^2 ^ + /34e2^h$2 

- Ei - B{p^ + B^Zl + E, 

+ ATTTe'y^iB^Zh + Ci - ^nTe^^Bpeij (j 


(2.34) 

(2.35) 


(2.36) 


where = ^{4’i'>’h))- Finally, the conductivities are obtained by differentiating the bound¬ 
ary currents (J*(oo), Q^{oo)) with respect to the external electric field(Fii) or the thermal 
gradient (0): 


^ dJ\^) ^ ^ + Z^ 6^^ , 


dEj 


dEi 


Sh) 


^B. 


1 dE{oo) , ^ if, ^ , 1 d5g)u 


1 dQ^ipo) 


M 


[h) 


= -47r5*^^ + , 


T dE, 

1 9Q*(oo) 


dE. 

M 


j 

[h) 


■4?' 

T' 




T dL 


dCj 


•4?' 

T‘ 


(2.37) 

(2.38) 

(2.39) 

(2.40) 


where we put hats on the conductivities to distinguish them from the ones where magneti¬ 
zation current are taken out [29]. More explicitly, with (2.36), the general DC conductivity 
formulas are given as follows. 

(i) The electric conductivity and Hall conductivity: 


(7 — — r, , 

(2.41) 

+ {B^Zh + 

-_X. -_.X Bp (p-^ + B^zl + 2e-^‘li^ZA) 

(2.42) 

B^p'^ + {B^Zh + e^'^/32$;,)2 
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(ii) The thermoelectric conductivity: 


cr 


XX _ ^yy _ 




a^y = -ay^ = 




+ {B^Zh + e^hi32<^^)^ 

Aire^’^B (p^ + B^Z^ + e^^f3^Zh^h) B 


B2p2 + + e^hp2^^Y 


+ -Si 

T 


(2.43) 

(2.44) 

(2.45) 


(iii) The thermal conductivity: 


^xx iiyy 
K — ^yy — 


-xy _ _gyx _ 


leTT^Te^’"* [B^Zh + 

+ {B-^Zh + e^'^/32$i,)2 
167r^Te^^^Bp 


B2p2 + ( 52 ^^ + e’'^; 02 ^^) 


2 +^^2 


(2.46) 

(2.47) 


where 


and 


S 2 


poo 

El = / (ir'e-’^(^')Z((/)(r')), 

Jf'h 

poo 

= —2 dr' a{r')e~'^^^'^ Z {(j){r')). 

yxh 


(2.48) 

(2.49) 


The yEj terms in (2.44) and (2.47) come from the contributions of the magnetization cur¬ 
rent and the energy magnetization current, which should be subtracted [29]. In particular, 
in the case of the dyonic black hole in section 3.2^, Si = ^ and S 2 = 2(Jvt ^ ^ 

is the magnetization and is the energy magnetization. The relation between Sj and 
the magnetizations are confirmed by (3.20). Therefore, the DC conductivities((T®'^, a*-^,/t*2) 
read 


, (2.50) 

d*^-|sie*^ (2.51) 

k*^-|s 26 *^ (2.52) 

which are expressed in terms of the black hole horizon data. 

2.1 Nernst effect 

The thermoelectric conductivites play an important role in understanding high Tc super¬ 
conductors. In the presence of a magnetic field, a transverse electric field can be generated 
by a transverse or longitudinal thermal gradient. The former is called the ‘Seebeck’ effect 
and the latter is called the ‘Nernst’ effect. 

The electric current, J, can be written in terms of the external electric field and the 
®The magnetization and energy magnetization current in a general setup were computed in [39] 
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thermal gradient as follows; 


(2.53) 


J = aE — aVT, 

where a and o. are 2x2 matrices. In the absence of the electric current, 

Ei = (2.54) 

Based on the definition of the Nernst effect, the Nernst signal (ew) is defined as 

eN =■ a)y^. (2.55) 


The Nernst signal in cuprates shows different features from conventional metals, so it is 
one of the important observables in understanding high superconductors. For example, 
in conventional metals the Nernst signal is linear in while in the normal state of a 
cuprate it is bell-shaped as a function of B. See, for example. Figure 12 in [35]. At a fixed 
B, the Nernst signal increases as temperature decreases in the normal state of a cuprate 
and, in turn, near the superconducting phase transition the Nernst signal becomes much 
stronger than conventional metals as shown in Figure 20 in [35]. 

Now we have the general formulas for the DC transport coefficients, (2.50) and (2.51), 
we can compute a general Nernst signal (2.55) 


pi + + Zl {B^p^ + ’ 


(2.56) 


which is expressed in terms of the black hole horizon data. By playing with the parameters, 
Zh,Vh, and <h/i, we can simulate ew by the formula (2.56). It may guide us in constructing 
more realistic models showing aforementioned cuprate-like properties and furthermore in 
understanding the physics of strongly correlated systems. 

There are two comments on general features of the Nernst signal (2.56). First, it is 
proportional to for small value of fi and it goes as 1//?^ for large value of /3. It has the 
maximum at /3 = /3max 


/3 


2 

max 


+ B'^Zl 
Zh^h 


(2.57) 


Second, in the limit, p —?• 0, the Nernst signal becomes 


cn 


p=0 


AttB 


(2.58) 


In this regime, relevant to the quantum critical point, the Nernst signal is proportional to 
the inverse of at fixed B and [3. Because the Nernst signal increases as temperature 
decreases in the normal phase of cuprates [35], to have a cuprate-like property, should 
decrease as temperature decreases. This will be a restriction on in model building. 
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3 Example: dyonic black branes with momentum relaxation 
3.1 Model with massless axions 

As an explicit model, we consider the Einstein-Maxwell system with massless axions. The 
action is given by ( 2 . 1 ) with the following choices 

0 = 0, $i = $2 = l, E(<^) = -6 /l2, Z(0) = 1, (3.1) 


where L is the AdS radius which will be set to be 1 from now on. Adding the Gibbons- 
Hawking term, we start with 


5o 


1 


IGttG 


JM 


R + Q- 





1 



d^x^/^K, 


(3.2) 


where 7 is the determinant of the induced metric on the boundary. K is the trace of 
the extrinsic curvature tensor Kmn‘^ defined by — 7 ^ 7 ^V(-pnQ), where n is the outward¬ 
pointing normal vector^ to the boundary(dM) which is at r = A. From here, for simplicity, 
we also take IGvrG = 1. For the holographic renormalization, we have to add a counter 
action 


s,= j dx^ ( -4 - i2[7] + J 

JdM \ ^ ^ 


•XI 


and the finite renormalized on-shell action is 


(3.3) 


Sr 


lim (S'o -b Sc) 

A^OO 


on—shell ‘ 


(3.4) 


Since the boundary terms do not change the equations of motion, the equations (2.2)-(2.5) 
are valid and yield, with (3.1), 


Rmn = ^QMN (r + Q- - - ^(< 9 x/)^ j + X ^ dMxANXi + i:Fm^Fnp , ( 3 . 5 ) 


1=1 


= 0 , 


vA/ = o. 


(3.6) 


We want to find a solution of the equations of motion, describing a system at finite 
chemical potential(^) and temperature(T) in an external magnetic field(i?) with momen¬ 
tum relaxation. It turns out the dyonic black brane solution modified by the axion hair 

^Where M,N — 0, l,2,r are the indices for bulk and fijU = 0,1,2 are the indices for the boundary 
coordinates. 

®In our case, = ^0,0,0, 1/y/U{r)^ . See (3.7). 
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(2.7) does the job. i.e. 


ds2 


= - + iTM 

U{r) = r'‘ - 


dr^ _l_ ^2 ^^^2 _|_ ^ 

mo + ql, r\ 

r ^ 4 J.2 ' 


a{r) = , B = qmrh , 

Xi = I3x, X 2 = I3y, 
where r/j is the location of the horizon and 


mo = rl 


f. I B^ + Qm 

1 H 



(3.7) 


(3.8) 


3.2 Thermodynamics 

To obtain a thermodynamic potential for this black brane solution, we compute the on-shell 
Euclidean action(S'^) by analytically continuing to the Euclidean time(r) of which period 
is the inverse temperature 


t = —IT , = —iSr 


(3.9) 


where S'® is the Euclidean action. By a regularity condition at the black brane horizon the 
temperature of the system is given by the Hawking temperature, 


T = 


U'jrh) 

dvr 




4r/i 


(3.10) 


and the entropy density is given by the area of the horizon. 


s = dvrr^ . 


(3.11) 


Plugging the solution (3.7) into the Euclidean renormalised action (3.4), we have the 
thermodynamic potential(H) and its density(W): 


qE '^‘2, o2 I 2 \ _ 

S = —(-mo-/3 rh + g^rh) = —W=- 


(3.12) 


where V 2 = / dxdy. The potential density W can be expressed in terms of the thermody¬ 
namic variables as 




W = -rl-^ + 2/3^ - 3g^) 


V 2 ■" 4 

= s — sT — HP ^ 


(3.13) 


where e and p are the energy density and the charge density respectively. The second line 
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is obtained by using the relation 


e = 2mo 


P = 


(3.14) 


which is derived as follows. We want to compute one-point functions of the boundary 
energy-momentum tensor(T^j^), current(J^) and scalar operators(0/) dual to xi- Our 
metric is of the following form 


-b x^iv{dx^ + V^dr){dx'^ + V'^dr) 


(3.15) 


where X and are the lapse function and the shift vector. The extrinsic curvature 
tensor has non-vanishing components, — D^Vy — DyV^j), where is a 

covariant derivative of the boundary metric Xtiv lu terms of aforementioned variables, we 
define ‘conjugate momenta’ of the helds as 


^ _ dSren 


= {K^y - Kx^y - 2Xf,y + Gf,y[x] - ^df,xidyXi + ' Vx/^ , 


= 


6Sren 

6Ay 


= , n/ = 


5Sr 


dxi 


= V=7(-w^x/- a^xi) 


(3.16) 


Thus, the one point functions are 


{Tf,y) = lim 


2r 


r^oc -»/ — 


= 


^ 2mo 0 0 

0 mo 0 
0 0 mo 


(J^) = lim 


r—^oo -v/ — 


n'" = ( prh, 0, 0 ) 


(0/) = lim 


r^oo A/ — 


= 0. 


(3.17) 


which yield (3.14). 

Since the pressure T = — W, (3.13) becomes a Smarr-like relation 

e + 7 = sT + . 

Notice that the pressure is not equal to (T^x) since 

7 = {Txx) + rh0^ - ruqli ■ 

The magnetization(M) and the energy magnetization(M^) are 


M = - 


f5n 


\6B 


— Qm : 


= - 


f 6X1 \ _ 1 


(3.18) 


(3.19) 


(3.20) 
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where SB and 6B^ are defined by a linearized solution® 


Sgty = xU{r)dBE , SAy = X i^B + ^SBe—'^ , 


B 


2r^ 


B 


SAt = {1-B^\5 Be - ^(1--]5B 




rh 


We find the first law of thermodynamics 

58 = TSs + gSp - rhS{p^) - MSB, (3.21) 

by combining the variation of the first line of (3.13) with respect to T,p,l3‘^ and B: 

5W = -MSB - sST - rhS{l3^) - pSp, (3.22) 


and the variation of (3.18). 


3.3 DC conductivities: Hall angle and Nernst effect 

In this section we study the DC conductivities of the dyonic black brane with momentum 
relaxation. Because we have derived the general formulas in section 2, we only need to 
plug model-dependent information (3.1) and (3.7) into (2.50)-(2.52). 

The electric conductivities yield 


^XX ^ ^yy _ p,2^2 




a 


xy _ 


= = Brhp- 


B"^ + rl + 2 / 5 ^) 


In the clean limit, fS 


a 


r\p?‘B‘^ -|- [B‘^ -|- 
0, these boil down to 

ThP P 


= ayy = Q, 


B 




(3.23) 

(3.24) 


where temperature dependence drops out and we recover the results expected on general 
grounds from Lorentz invariance, agreeing with [28]. In the limit i? —?■ 0, the expressions 
become 


„xx vy 1 I 
cr = = f + 7^ 


a 


xy _ 


= ay^ = o, 


®For details we refer to Appendix C of [29]. 
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Figure 1. /3/T and B/T dependence of DC conductivities at fixed ^/T = 4 


which reproduces the result in [16, 19]. The thermoelectric and thermal conductivities read 


^xx ^yy 
a = = 


4TTrll3‘^fi 


v2 ’ 


a^y = —ay^ = AirriB 


rln^B^+ {B‘^ + rll3‘^Y 


rlY^B^ + (52 + 


2 ’ 


= Kyy = lOvrrfT- 


B^ + r2/32 


R^y = -kV^. 


+ {B^ + rlP^) 

IGn'^rf^RiTB 

r2^252+ (52+^2^2)2 ■ 


2 ’ 


(3.25) 

(3.26) 

(3.27) 

(3.28) 


To see the effect of /3 and B on the conductivities, the formulas (3.23)-(3.28) are not 
so convenient since rh is a complicated function of T,B,rl, and (3, as expressed in (3.10). 
Therefore, we make plots of conductivities in Figure 1, where we scaled the variables by 
T and fixed rl/T = 4. The and are qualitatively similar; the B dependence 

at fixed /3 is monotonic, while the f3 dependence at fixed B is not. The a^y,a^y, and R^y 
are similar; the /3 dependence at fixed B is monotonic, while the B dependence at fixed 
(3 is not. See Figure 4 for the cross sections at B/T^ = 1. As the chemical potential 
increases, the electric and thermoelectric conductivity generally increases while its overall 
2-dimensional shape does not change qualitatively. However, the thermal conductivities 
behave differently; R^^ broadens while R^y is sharpened. 
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0 


2 




(a) B-P dependence at /i/T = 4 (b) /3-T dependence a.t B/fi — 1 


Figure 2. Hall angle 9h = arctan(cr^^/cr^“) 


(a) P/T and B/T^ dependence of 
the Nernst signal at fj,/T = 1 





(b) The cross sections of (a) at 
P/T = 0.5, l,4(red, green, blue) 


Figure 3. Nernst signal 



The Hall angle(0//) is defined by the ratio of the electric conductivities: 


isnOu = 


a' 


xy 




^jlB + 2 / 3 ^) 

+ r 2(^2 


(3.29) 


which agrees to the result reported in [23]. As shown in Figure 2, the Hall angle 9h ranges 
from 7r/2 to 0. The angle increases as B increases or (3 decreases(Figure 2(a)). In the 
strange metal phase, we are interested in the temperature dependence of the Hall angle, 
which is proportional to 1/T^. In our case, we numerically found that the Hall angle ranges 
between 1/T^ and l/T^(Figure 2(b)). In the large T regime, the Hall angle always scales 
as 1/T. It can be seen also from the formula (3.29), where if T is large compared to the 
other scales, r// ~ T so ianOn ~ l/T. 

The Nernst signal (2.56) yields 


^ 47rrfe/3 B 
^ ^ 2^2 +^ 2(^2 + ^ 2 ) 2 - 

As discussed in Section 2.1, the Nernst signal is linear in B in conventional metals while 
it becomes bell-like in the normal state of cuprates [35]. To see whether our model can 
capture this feature, we make a three dimensional plot of the Nernst signal as a function of 
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/3/T and BjT’^ at fixed ijl/T = 1 in Figure 3, where (b) is the cross section of (a) at fixed 
(5. The blue line is almost straight while the green and red ones are bell-like. Therefore, we 
find that our system shows the transition from the normal metal(blue line) to cuprate-like 
state(green and red) as j3 decreases, in the Nernst signal perspective. The green and red 
curves are similar to Figure 12 in [35] and it was proposed that the bell shape can be 
explained by the dynamics of the vortex liquid in non-superconducting phase [36, 37]. It 
was also interpreted as an evidence of the pseudo-gap phase [35]. If we take this point of 
view, our model might be relevant to the pseudo-gap phase. 

One may wonder if the blue line is qualitatively similar to the green and red ones 
for higher B/T‘^. As B/T'^ increases the blue line reaches the maximum and has almost 
plateaued out, which is not bell-shaped. Furthermore, the numerical value of B/T'^ at the 
maximum eAr(~ 60) is about 1600, which is too big to consider, compared to other scales. 

4 Numerical AC conductivities 

So far we have discussed the DC conductivities. In this section, we will consider the AC 
conductivities. To be concrete, we continue to investigate the model in the previous section, 
namely, the dyonic black brane with the momentum relaxation by the axion fields. 

4.1 Equations of motion and on-shell action 

In section 2 and 3, we have derived the DC conductivity in a gauge gri / 0. In this section, 
to compute the AC conductivity, we will work in a gauge 6gri = 0. The physical result 
should not be changed by this gauge choice^. Indeed, we will show the zero frequency limit 
of the AC conductivity (with 6gri = 0) agrees to the DC conductivity (with 6gri / 0) in 
section 2. It serves as a good consistency check of our numerical method. 

To compute conductivities holographically it is consistent to turn on linear fluctuations 
of bulk fields, dgu, SAi and bxi at zero momentum in the x and y directions, where 
i = I, 2(or i = x,y). The fields 6gu and SAi are related to the heat current and the electric 
current. They were introduced in [28] and here we add 5xi for momentum relaxation. The 
fields, 5gti, SAi and Sxi, can be expressed in momentum space as 

/ OO J 

—e—*a,(a;,r), (4.1) 

/ OO J 

(4.2) 

-OO 

/ OO J 

V.(a;,r), (4.3) 

where in the metric fluctuation (4.2) is introduced to make the asymptotic solution of 
hti constant at the boundary(r —>• oo), for the sake of convenience. From (3.5) and (3.6), 
the linearized equations for the Fourier components are given as follows. 

^There is a subtle issue on the gauge choice and holographic renormalization. For more details we refer 
to [33]. 


- 16 - 



- Einstein equations: 


ghhu iujqmaj P'^hu fia'i 4:h[^ „ _ 

r^U r'^U r'^U r'^U r i* ’ 


iUqma'j iixQmhtj 

- 4 - 1 - 2 -^ - 4 - 

r*oj r^U! 


+ ^ + Ki —0 ) 


- Maxwell equations: 


t/'a' iuqmhtj u'^ai „ 

U ^ [/2 + [/ + ^72 


- Scalar equations: 


17V' ipujhti 

U 


Among these eight equations, only six are independent. 

Near the black hole horizon (r —)■ 1)® the solutions are expanded as 


V = (r-ir±+HhSP+hf)(r-l) + ...), 

Oi = (r - + af^^(r - 1) + • • •), 

V'i = (r - 1)''± (^ - 1) H-) > 


(4.4) 

(4.5) 


(4.6) 


(4.7) 


(4.8) 


where u± = ±i4a;/(—12 + + 2,0^ + /i^). In order to impose the incoming boundary 

condition relevant to the retarded Green’s function [41], we have to take v = It turns 
out that the 4 parameters a^P and ipP may be chosen to be independent since hP and 
all higher power coefficients can be determined by them. 

Near the boundary(r —>■ oo), the asymptotic solutions read 


hti = hP + ^hP + ^hP + ■ ■ 
ai = aP + ^aP + • • • , 

+ • • 


(4.9) 


where the leading terms hP,aP, and are independent constants, which fix 

completely by the equations of motion. The nontrivial sub-leading terms, hP,aP, and 
f3) 

' will be determined by the incoming boundary conditions at the horizon for the given 


-,(o) 


.(2) ,/,(2) 


leading terms, hP,aP, and ipr^- leading terms play the role of sources for the 

operators of which expectation values are related to hP,aP, and pp respectively. 

Expanding the renormalized action (3.4) around the dyonic black brane background 


-,( 0 ) 


®From here we set = 1 for convenience of numerical analysis. 
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and using the equations of motion, we obtain a quadratic on-shell action: 



1 f o 

lim - / d X 
A^oo 2 Jr=A 


4r^ - 


4^4 




fhl - 


Pr „ „'r r" 


VUir 


-SxiShti 






-U{r)6Ai5A'i -h r^6hti5h[i - r‘^U{r)5xi^x!i Shi ( fjAAi - 


Pr'^Sxi 

VW) 


(4.10) 


where dhu = r~'^Sgu{t,r), ‘dot’ denotes dt and ‘prime’ denotes dr- We have dropped the 
contributions from the horizon, which is the prescription for the retarded Green’s function 
[41], In particular, with the spatially homogeneous ansatz (4.1)-(4.3), the quadratic action 
in momentum space yields 



V 2 (S 

2 Jo 27r 


"'ti 


2moh[fhlf + - 3h[fh[f + 3 


(4.11) 


where V 2 is the two dimensional spatial volume f dxdy. The argument of the variables 
with the bar is —u. We dropped the range (—oo,0), which is the complex conjugate of 
(4.11), to obtain complex two point functions [41]. 

The on-shell action (4.11) is nothing but the generating functional for the two-point 
Green’s functions sourced by We may simply read off a part of the two 

point functions from the first two terms in (4.11). The other three terms are nontrivial 
and we need to know the dependence of {a\^\ However, the 

linearity of the equations (4.4)-(4.7) makes it easy to hnd out a linear relation between 
and {af'\ In the following subsection we will explain how to 

find such a relationship numerically in a more general setup and apply it to our case. 


4.2 Numerical method 

A systematic numerical method with multi fields and constraints was developed in [19] 
based on [42, 43]. We summarize it briefly for the present case and refer to [19, 33] for 
more details. Let us start with N fields 4>“(x, r), (a = 1, 2, • • • , N), which are fluctuations 
around a background. Suppose that they satisfy a set of coupled N second order differential 
equations and the fluctuation fields depend on only t and r: 


4>“(t,r) 


dw 

W) 


e-^-hPK(r), 


(4.12) 


where is multiplied such that the solution 4>“(r) goes to constant at the boundary(r —>■ 
00 ). For example, p = 2 in (4.2). 

Near horizon(r = 1), solutions can be expanded as 


4>“(r) = (r - 1)''“± ((/?“ -h 0“(r - 1) H-) , 


(4.13) 


where we omitted the subscript w for simplicity and Ua+{i^a-) corresponds to the incom- 
ing(outgoing) boundary condition. In order to compute the retarded Green’s function we 
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choose the incoming boundary condition [41], This choice reduces the number of indepen¬ 
dent parameter from 2N to N. There may be further reductions by Nc if there are Nc 
constraint equations. As a result, the number of independent parameter is — A^c so we 
may choose N — Nc initial conditions, denoted by = 1,2,-- - , N — Nc): 


¥^“3 


^ N-Nr 


(\ 1 1 ...\ 
1-1 1 ... 

1 1 - 1 ... 

Vi 1 1 ...) 


(4.14) 


Every column vector yields a solution with the incoming boundary condition, denoted 
by which is expanded as 

O? 

> S? -h • • • H— j- + • • • (near boundary), (4.15) 


where are related to the sources, which are the leading terms of Tth solution, and 
are related to the expectation values of the operators corresponding to the sources(5a > 1). 
A general solution constructed by <h“j(r) is 


- - 0?c* 

4)? (r) = <l)?(r)c* S?c* 4-^-^-^- (near boundary), (4.16) 

2 2 


with real constants c*. We want to identify S?c* with the independent sources J“ but if 
there are constraints, it is not possible since a > i. However, in this case, there may be Nc 
other solutions corresponding to some residual gauge transformations [19, 33] 

0 ? 

^ Sf + • • • H —^ + ■ ■ ■ (near boundary), (4.17) 


where i runs from A^ — A"c + 1 to A^. This extra basis set, <h“(r) = 4?“c*, generates a general 
solution together with the ingoing solutions. 

In our case, A^ = 6 and Nc = 2, which corresponds to the constraints pri = 0. There 
are two sets of additional constant solutions of the equations of motion (4.4)-(4.7) 


hfi — h 


ti ) 


iqmeijh^tj 


Xi = 


il3hl 


ti 


(Li — 




UJ 


(4.18) 







where is arbitrary constant and i,j = 1,2. Therefore, the explicit expression for is 


S“.= (S“5 §“6 


/ 0 

UJ 


: Qm 

■' U) 

1 




\ 0 


0 
0 
1 
0 

ip ) 

CaJ / 


These can be understood as residual gauge transformations keeping gri = 0, which are 
generated by the vector fields, of which non-vanishing components are and 

constants), i.e. L^Ai = ^^9ti = and = 

Therefore, the most general solution reads 

^rn(^) + ^c(^) = + - + ^ + (4-19) 

where we defined J“ and For arbitrary sources J“ we can always find 


c 


I 




(4.20) 


where I = 1,... ,N. The corresponding response is expressed as 

R'^ =0^jc^ 


(4.21) 


The general on-shell quadratic action in terms of the sources and the responses can be 
written as 



Va r°^ da; 

^ Jo W) 


J^KbJ^ + nabR^ , 


(4.22) 


where A and B are regular matrices of order N and the argument of J“ is —oj. For example, 
the action (4.11) is the case with: 








ay 


ay 

'hx 

'^ty 

, = 

'hx 

'^ty 

i>f 


A' 

At"/ 


At'/ 


A = 

/O -g 0\ 
0 —2mo 0 

® l2, 

B = 

A 0 o\ 

0-3 0 


o 

o 

o 



o 

o 


(4.23) 

where the index u) is suppressed and la is the 2x2 unit matrix. Plugging the relation 
(4.21) into the action (4.22) we have 




dw 

W) 




(4.24) 
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which yields the retarded Green’s function 


Gab — ^ab ■ 


(4.25) 


In summary, to compute the retarded Green’s function, we need four square N x N 
matrices. A, B, § and O. The matrices A and B can be read off from the action (4.22), which 
is given by the on-shell expansion near the boundary. The matrices S and O are obtained 
by solving a set of the differential equations. Part of them comes from the solutions with 
incoming boundary conditions and the others may be related to the constraints. Notice 
that the Green’s functions do not depend on the choice of initial conditions (4.14). 

In the case of the dyonic black branes in section 4.1, we may construct a 6 x 6 matrix 
of the retarded Green’s function. We will focus on the 4x4 submatrix corresponding to 
and in (4.9), 


gIj 

^TJ 


^.JT 


(4.26) 


where every is a 2 x 2 retarded Green’s function with i = x,y for given a and (5. The 
sub-induces a, (3 denote the operators corresponding to the sources, i.e. The is dual 
to the electric current J* and the is dual to the energy-momentum tensor T**. From 
the linear response theory, we have the following relation between the response functions 
and the sources: 


( ("'•> ) ^ ( a‘jj a% \ f af > \ 

({T«)j [cx’j ) [hlf ) 


(4.27) 


where (J*), (T**),a^-°^ and /i® are understood as 2 x 1 column matrices, with i = x,y. We 
want to relate these Green’s functions to the electric((T), thermal(«:), and thermoelectric(Q;, a) 
conductivities defined as 


( {.r)\ ^ a^^T\ ( E, \ 

yiQ^) J ) y-{VjT)/T ) ' 


(4.28) 


where Q* is a heat current, Ei is an electric field and VjT is a temperature gradient along 
the i direction. Notice that the electric and heat current here contain the contribution of 
magnetization, so we use the conductivities with hat((2.37)-(2.40)). By taking into account 
diffeomorphism invariance [1, 2, 19], (4.28) can be expressed as 


( (J*) \ 63^t\ (iuj{af + yLhfh\ 

[ (T^*) - y{E) j “ ) [ iu;hlf ) ' 


(4.29) 


From (4.27) and (4.29) with the magnetization subtraction (2.50)-(2.52), the conductivities 
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Re[<T“] Re[a“l Re[it“] 




Figure 4. Agreement of DC analytic formulas(solid curves) obtained in section 3.3 with numerical 
results(red dots): DC conductivities vs ^/T with /x/T = 4, B/T"^ = 1 


are expressed in terms of the retarded Green’s functions as follows 




iGlr 


G^j) 



B f 0 

T j ' 

(4.30) 


4.3 AC conductivities and the cyclotron poles 

In this section we present our numerical results. Some examples of the AC electric con¬ 
ductivity are shown in Figure 5 and 8; the thermoelectric conductivity is in Figure 6; and 
the thermal conductivity is in Figure 7. 

Before discussing the AC nature of the conductivities, we first examine the DC limit(w —)> 
0) of the AC conductivities and compare them with the analytic expressions derived in sec¬ 
tion 3.3. The comparisons are shown in Figure 4, where all conductivities are plotted as a 
function oi j3/T with the other parameters fixed; /x/T = 4 and B/T^ = 1. The solid lines 
were drawn by the analytic expressions, (3.23), (3.24), (3.25)-(3.28) and the red dots were 
read off from the AC numerical results in the limit a; —)> 0. Both results agree, which serves 
as a supporting evidence for the validity of our analytic and numerical methods. Indeed 
this agreement is not so trivial in technical perspective. In the DC computation we turned 
on hri and read off the physics from the horizon data while in the AC computation we 
worked in the gauge hri = 0 and considered full evolution in the r-direction. 

Let us turn to the AC properties of the conductivities. Figure 5, 6, and 7 show 
the /? dependence of the electric, thermoelectric, and thermal conductivity respectively. 
The dotted lines are the cases for /? = 0 and the red, green, blue curves are for fi/T = 
0.5,1,1.5. The red dots at w = 0 in the real part of the conductivities are the values 
given by the analytic formulas (3.23), (3.24), (3.25)-(3.28). They agree to the numerical 
AC conductivities in the limit w —>• 0. As /3 increases, when all other scales are hxed, the 
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(b) Im ( 


(c) Re 


0.1 0.2 0.3 0.4 0.5 

(d) Im 


Figure 5. /3 dependence of electric conductivities: /3/T = 0, 0.5, 1, 1.5(dotted, red, green, blue) 
with ixjT = 1, 5/T2 = 3 
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(a) Re 


(b) Im 


ReEa"’] 



o.i'. <31 .0:3 0.4 0.5 



(c) Re 


0.1 0.2 0.3 0.4 0.5 

(d) Im 


Figures. /3 dependence of thermoelectric conductivities: /3/T = 0, 0.5, 1, 1.5 (dotted, red, green, 
blue) with fi/T — 1, B/T"^ = 3 




0.1 0.2 0.3 0.4 0.5 
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(a) Re k"" 


(b) Im 



(c) Re 


0.1 0.2 0.3 0.4 0.5 

(d) Im 


Figure 7. /3 dependence of thermal conductivities: /3/T = 0, 0.5, 1, 1.5(dotted, red, green, blue) 
with ^l|T= 1, B/T'^ = 3 
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(c) Re 


0.1 0.2 0.3 0.4 0.5 

(d) Im 


Figure 8. B dependence of electric conductivities: B/T^ = 0.5, 1, 2, 4(red,orange,green,blue) 
with j3/T = 1/2, fi/T = 4 


curves become flatter, which is expected from stronger momentum relaxation. Notice that 
there is no l/cj pole in the imaginary part of the condnctivities even when /? = 0. It is 
becanse the gauge field for B (2.6) breaks translation invariance in the same way as the 
axion fields do. 

There is a peak in the cnrves in Figures 5-7, which is related to the cyclotron resonance 
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pole, the position of the pole of the conductivity in complex oj plane [29, 30] 


OJ^: = UJc — , 


(4.31) 


where the cyclotron frequency ojc is the relativistic hydrodynamic analog of the free particle 
case, w/ = eBjme. However, the resonance due to w* here should be understood to be due 
to a collective fluid motion rather than to free particles. A damping 7 could be thought 
of as arising from interactions between the positively charged current and the negatively 
charged current of the fluid, which are counter-circulating. In the hydrodynamic regime 
for small H at /? = 0, the w* was computed analytically in [30] as follows: 


_ 0 ^ — 0 


(4.32) 


Because j3 is related to momentum relaxation, it is expected that 7 will increase as f3 
increases. It is indirectly shown in the plots since all curves become flat, which may reflect 
the fact the pole goes away from the real u axis. It turns out that the ujc tends to increase 
as /3 increase, although it is not so clear in the plots. It will be discussed later in Figure 10 
and the equation (4.34). While in Figure 5-7 we focused on the effect of f3 at fixed B, in 
Figure 8 we investigated the effect of B at fixed (3; the red, orange, green, and blue curves 
are for B/T'^ = 0.5,1, 2,4 respectively. A relatively big /i//3 = 8 has been chosen since the 
peak is shaper when n/(3 is bigger. We find the peaks of the curves shift towards higher 
frequencies as B increases, which is consistent with the hydrodynamic analysis (4.32). 

To see the cyclotron pole directly in complex uj plane, we introduce the following 
combination 

± , (4.33) 

for easy comparison with [29, 30]. The density plots of jci+j in complex u plane are shown 
in Figure 9. Figure 9(a),(b),(c) are the cases with j3 = 0 and 9(d),(e),(f) are the cases with 
(3/T = 10. We choose fj. and Qm such that to be constant. First, at /3 = 0 we 

recover the result of [30]®. White areas correspond to poles and dark areas are zeroes of . 
There is a symmetry: fj. ^ and q^ —t at /3 = 0 inherited from the electromagnetic 

duality of the bulk theory [30]. Since this bulk duality holds at finite j3 we expect the 
same symmetry is preserved. It is demonstrated in Figure 9(d) (e)(f); the figure (d) and 
(f) are symmetric under the exchange of the white and dark region. The finite (3 shifts the 
position of the poles to the negative imaginary direction. This implies the width of the 
peak increases at real uj axis as discussed previously. 

The magnetic field dependence of the cyclotron poles at different values of f3 is shown 
in Figure 10, where (a) and (b) show the real part(wc) and the imaginary part(— 7 ) re¬ 
spectively. The gray, red, green, and blue dots are the numerical results for /3 = 0, 2, 3, 4 
respectively, while the black dashed line is the analytic result at /3 = 0 for small B from 
the hydrodynamic analysis (4.32). At /3 = 0, our numerical result(gray dots) agrees to [30] 
and also fits well to the hydrodynamic analysis (black dashed line) for small magnetic field. 

®Our numerical values of fi and qm are different from [30] due to convention difference. 
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(a) qmjT = 0, iijT = 12.5 


(b) q^jT = 8.9, iijT = 8.9 (c) q^jT = 12.5, iijT = 0 




(d) q^/T = 0, ti/T = 12.5 


(e) qm/T = 8.9, ^/T = 8.9 


(f) q™/T = 12.5, M/r = 0 


Figure 9. A density plot of |cr"''| in complex w plane: White areas correspond to poles and dark 
areas are zeroes of ct"*". /3 = 0 for (a)(b)(c) and (i/T = 10 for (d)(e)(f). 


Re[ai.] 


(a) Re uj,{=uJc) 


lm[b;,] 



(b) Im ujt{= — 7 ) 



Figure 10. The magnetic field dependence of the position of the cyclotron resonance pole: /3 = 
0, 2,3,4(gray, red, green blue). The dotted lines are the results by the hydrodynamic analysis (4.32). 


We investigated how w* changes as f3 is introduced. Based on the red, green, and blue dots 
in Figure 10 and additional similar numerical data for different parameters, we found the 
following relation at small B 

= Uc - ij = uJc + ci0^B - f(7° + C2/3^), (4.34) 

where and 7 ® are defined in (4.32). It seems that ci and C 2 are dimensionful constants 
independent of (3 and B. For this formula we focused on small B region where the dots in 
Figure 10(a) are linear to B. However, at large magnetic field B we numerically found the 
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tendency that ujc ~ csi?, where C 3 seems independent of fi. 

In the presence of dissipation, the cyclotron frequency was shown [29] to be changed 
as 

7^7 + , (4.35) 

'^imp 

while LJc is intact. In our case Timp is proportional to l//3^ for small j3 so the shift of the 
imaginary part in (4.34) is consistent with the hydrodynamic calculation [29]. However, 
our result implies that the cyclotron frequency(wc) is also shifted by ciP^B ~ B/rimp- We 
suspect that the analysis in [29], where B/T^ <C 1 is assumed, is valid in the limit ciB is 
small. We leave this issue and the analytic justification of the specific form (4.34) for a 
future project. 

5 Conclusions 

In this paper, we have computed the electric, thermoelectric, and thermal conductivity 
at finite magnetic field by means of the gauge/gravity duality. First, by considering a 
general class of Einstein-Maxwell-Dilaton theories with axion fields imposing momentum 
relaxation, we have derived the analytic DC conductivities, which are expressed in terms 
of the black hole horizon data. As an explicit model we have studied the dyonic black 
hole modified by a momentum relaxation effect. The background solution is analytically 
obtained and the AC electric, thermoelectric, and thermal conductivity were numerically 
computed. The zero frequency limit of the numerical AC conductivities agree to the DC 
formulas. This is a non-trivial consistency check of our analytic and numerical methods 
to compute conductivities. Our numerical method can be applied to other cases in which 
multiple transport coefficients need to be computed at the same time. 

The Nernst signal, the Hall angle, and the cyclotron resonance pole were discussed 
following [28-30]. Our general analytic formulas of the Nernst signal can be used to build 
a realistic model and to investigate the universal properties of the model. In particular, 
in the dyonic black hole case, the Nernst signal is a bell-shaped curve as a function of the 
magnetic field, if momentum relaxation is small. It is similar to an experimental result 
in the normal state of cuprates [35], which was speculatively explained by a vortex-liquid 
effect [36, 37]. For large momentum relaxation the Nernst signal is proportional to the 
magnetic field, which is a typical property of conventional metals. The Hall angle for the 
dyonic black hole was computed explicitly. The Hall angle ranges between 1/T° and 1/T^ 
and scales as l/T for large T. However, in the strange metal phase, it was known that the 
Hall angle is proportional to 1/T^. The cyclotron poles (w*) we found are consistent with the 
hydrodynamic results at /3 = 0 [29, 30]. They are shifted by momentum relaxation(/3 7 ^ 0) 
and our numerical analysis suggests the specific dependence of the cyclotron poles on B 
and /I, (4.34), when B is small. If B is large it was proposed that cOc ~ C 3 B with C 3 
independent of j3. We plan to investigate the properties of the cyclotron poles in more 
detail both numerically and analytically. 

It is important to compare our AC conductivities with the general expressions based 
on the memory matrix formalism [40, 44]. It would be also interesting to compare our AC 
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conductivity results with [45], where the AC electric conductivities have been studied at 
finite magnetic field in the probe brane set up, focusing on the transport at quantum Hall 
critical points [46]. The metal phase of our model at H = 0 does not have the property of 
linear-T resistivity, so it is worthwhile to start with the models having that property and 
then investigate the Hall angle and the Nernst effect in those models. 
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